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Abstract 

We present a method to construct random model power grids that closely match 
statistical properties of a real power grid. The model grids are more difficult to 
partition than a real grid. 

1 Introduction 

Power grids are prime examples of the interconnected networks that make life 
in technological societies possible. It is therefore of paramount importance 
to develop methods to prevent cascading failures that manifest themselves in 
widespread, catastrophic blackouts. One such defensive strategy is Intentional 
Intelligent Islanding [1] [2] , which aims to prepare for the rapid isolation of parts 
of the grid where instabilities arise before they can spread further. This problem 
of network partitioning [3, 4 involves determining subdivisions of the grid that 
are tightly internally connected, but weakly connected to the rest of the grid. 
At the same time, each such island should be close to self-sufficient with power. 

In previous papers [5j [6] we have considered network-theoretical methods to 
achieve this goal, using the high-voltage network in the U.S. state of Florida as 
a test example. Here we continue this endeavor by constructing model power 
grids that share important statistical properties with a real one. Such models 
provide opportunities to study the effects of specific network modifications, as 
well as to perform scaling analyses in terms of the network size. 
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2 Model construction 



The Florida high-voltage power grid [7] consists of N = 84 vertices, 31 of which 
are generators and the rest are distribution substations that act as loads. The 
vertices are connected by M = 200 edges, some of which are parallel power lines 
connecting the same two vertices. We use a simplified representation of the grid 
as a weighted, undirected graph [3], defined by the N x N symmetric weight 
matrix W, whose elements wtj > represent the "conductances" of the edges 
(transmission lines) between vertices (generators or loads) i and j, 

number of lines between vertices i and j 
lJ normalized geographical distance ' 

where the "geographical distance" is the length of the edge connecting i and 
j. To obtain a model independent of any specific systems of length units, the 
distances are normalized such that the areal density of vertices is unity. In 
Fig. Q] we show a map of the Florida network together with a representative 
model network. 

Model networks were produced by the following procedure. 

1. We placed the N = 84 vertices randomly in a square of area N. 

2. Following the standard "stub" method jS], we attached 2M = 400 stubs 
or half-edges randomly to the N vertices. (Actually, to ensure that no 
vertices in this small network should be totally isolated, we first attached 
one stub to each vertex, and then distributed the remaining 2M — N stubs 
randomly between the vertices.) The resulting degree distribution for the 
particular model grid discussed in this paper is shown together with that 
of the real Florida grid in Fig. [2{a). 

3. We connected the stubs randomly in pairs, with the restriction that self- 
loops (two mutually connected stubs at the same vertex) were forbidden. 

4. To obtain an edge-length distribution with the same average as that 
of the real Florida grid (« 1.09 in our dimensionless units), we em- 
ployed a Monte Carlo (MC) "cooling" procedure using a "Hamiltonian" 
in which the total edge length L plays the role of the system energy. 
The update attempts consisted in choosing two different edges, ij and kl 
with i 7^ j 7^ k =/= I, interchanging j and I, and calculating the change 
AL = [L(ij) + L(kl)\ — [L{il) + L(kj)\. Attempts were accepted with 
the Metropolis probability with a fictitious "temperature" T, P(ij, kl — > 
il, kj) = Min[l, exp(— AL/T)]. In the limit that all edges are <C y/N, it is 
easy to show that the partition function for this model is Z — 2ttT 2 , and 
so (L(ij)) = 2T. Edge length distributions before and after "cooling" are 
shown together with that of the real Florida grid in Fig.[2][b). The average 
edge length vs. the number of MC steps is shown in Fig. 02 together with 
a schematic of the MC update mechanism. 
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3 Network partitioning 



The method used to partition the grids was described in Ref. [6] and is only 
sketched here. 

Agglomeration. A trial partitioning is obtained by associating each load i with 
its "nearest" generator j, defined as the one to which the effective resistance 
[3], Rjj is minimum. 

Optimization. The goal is to obtain a partitioning into K islands, C = 
{Ci, Ck }, that balances the requirements for islands that are (i) strongly 
connected internally, but sparsely connected to each other, and (ii) approxi- 
mately self-sufficient with power. 

Success in achieving requirement (i) is measured by the modularity Q [3] . It 
compares the proportion of edges internal to islands with the same proportion 
in an average null-model. 

w V w J 

v 

where Wi = J2j w ij> w = Tlij w iji an d & C(j)) — 1 if vertices % and j 

belong to the same island, and otherwise. (Other quality measures could also 
be used [TO].) 

We show in Ref. [6] that success in achieving requirement (ii) requires min- 
imizing the sum of the squares of the currents Ii that enter or leave the grid 
at the individual vertices (positive for generators and negative for loads). As 
shown in Fig. 01 the generating capacities of power stations are highly correlated 
with their degrees, fcj. Assuming that the whole grid is in power balance, we 
approximate Ii = ^ — r- for generators and Ii — — ^ ki , for loads. 

^generators K J Z^loads K 3 

Weighting the two requirements equally (other weighting choices could be 
made) , we attempt to maximize by MC simulated annealing the quality measure 



Qinit V Ei-^ 2 )init) 



where the subscript "init" designates the value after the first recombination, but 
before any MC steps. The MC steps consist in moving loads that are peripheral 
to one island to a neighboring island to which it is also connected. 
Iteration. The optimized islands form a new network (analogous to real-space 
renormalization-group calculations), in which each island is represented by a 
vertex. The connections between the new vertices are the same as those between 
the previous islands. This defines a new conductivity matrix. The components 
of the new current vector, \I) represent the generating surplus or deficiency of 
each of the old islands or new "super-generators" (smooth symbols in Fig. [TJ 
or "super-loads" ("geared" symbols in Fig. [1]), respectively. This process of 
agglomeration and optimization is iterated until all the original vertices belong 
to one island, and the optimum partitioning is identified. The quality measure 
E is shown vs. accepted MC step in Fig. 03a) for the real Florida grid, and for 
our representative model system in Fig. EJb). 
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4 Results and Conclusion 



The model power grids introduced here were constructed to match the size, 
proportion of generators, average degree, and scaled edge length of the real 
Florida high-voltage grid. In fact, both the degree distribution and the full 
edge-length distribution for the two networks are quite similar. Nevertheless, it 
is easier for our partitioning algorithm to find a partition with a high value of E 
for the real Florida grid than for the models. The resulting partitionings, shown 
in Fig. [1] consist of eight islands with E max s» 0.8 for Florida and 17 islands 
with E max « 0.4 for the representative model. For Florida, the four largest 
islands comprise 75 of the 84 vertices. The "North- West" portion of the model 
grid appears particularly difficult to partition. We believe this indicates that 
the real power grid is more strongly correlated than the randomized models. 
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Figure 1: The Florida high- voltage power grid (top) and a representative 
model network (bottom). Generators are represented by squares and loads by 
ovals. The partitions shown are in each case the "best" ones obtained in ten 
independent runs with the bottom- up partitioning algorithm of Ref. [5] . Islands 
arc identified by different colors and by numbers that can be seen if viewed at 
high magnification. See text for discussion. 
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Figure 2: (a) Degree distributions for the Florida power grid (left set of bars, 
in red) and the representative model grid analyzed here (right set of bars, in 
blue). Also shown is a Poisson distribution for k — 1 with mean 2M/N — 1. (b) 
Cumulative edge-length distributions, normalized to M = 200. Ten independent 
realizations of the model are included (thin curves in the background), together 
with their average distributions (heavy curves in the foreground). The right- 
hand set of curves represents the initial, "infinite-temperature" distributions 
(their average gives an estimate for the model's density of states). The left- 
hand set of curves represents the final, "low-temperature" distributions at T — 
1.09/2 = 0.545. The bold, blue curve on the left represents the distribution for 
the real Florida grid. Its average edge length is approximately 1.09. The dashed 
parabola represents the theoretical model density of states for an infinite system 
with unit areal vertex density and average degree 2M/N. It fits the simulation 
well for short edge lengths. 
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Figure 3: Average edge length vs. MC steps per edge. Ten independent sim- 
ulation runs (thin curves in the background) and their average (thick curve in 
the foreground). Inset: schematic of the MC update ij, kl (blue) — > il, kj (red). 
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Figure 4: Generating capacities of Florida power plants vs. their degrees. The 
line is a least-squares fit to the data. 
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Figure 5: Quality measure E vs. accepted MC steps in the partitioning algo- 
rithm for Florida (a) and our representative model system (b). The black dots 
represent the maximum value of E obtained in each stage of the partitioning 
algorithm. The global maxima in both cases represent the highest E max ob- 
tained over ten independent partitioning runs. These maxima correspond to 
the partitionings shown in Fig. [T] 
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